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We have carried out a comparative study of the electronic specific heat and electronic structure of 
a and 5-plutonium using dynamical mean field theory (DMFT). We use the perturbative T-matrix 
and fluctuating exchange (T-matrix FLEX) as a quantum impurity solver. We considered two 
different physical pictures of plutonium. In the first, 5/ J + , the perturbative treatment of electronic 
correlations has been carried out around the non-magnetic (LDA) Hamiltonian, which results in an 
f occupation around a bit above n/ = 5. In the second, 5/ 6 — , plutonium is viewed as being close to 
an 5/ 6 configuration, and perturbation theory is carried out around the (LDA+U) starting point 
bit below nf =6. In the latter case the electronic specific heat coefficient 7 attains a smaller value 
in 7-Pu than in a-Pu, in contradiction to experiment, while in the former case our calculations 
reproduce the experimentally observed large increase of 7 in <5-Pu as compared to the a phase. This 
enhancement of the electronic specific heat coefficient in <5-Pu is due to strong electronic correlations 
present in this phase, which cause a substantial increase of the electronic effective mass, and high 
density of states at Ef- The densities of states of a and 5-plutonium obtained starting from the 
open-shell configuration are also in good agreement with the experimental photoemission spectra. 

PACS numbers: 



I. INTRODUCTION 

Plutonium metal is arguably the most complex 
elemental solid. It exhibits a large number of al- 
lotropes, and in all the phases its Pauli-like mag- 
netic susceptibility and resistivity are an order of 
magnitude larger than for simple metals. The phase 
transitions between different phases are sometimes 
accompanied by very large discontinuities in volume 
[l|. There is a generalized suspicion that the origin 
of the unique properties of Pu must be connected to 
its unique position in the actinide series lying at the 
boundary between the early actinides where the f 
electrons are band-like and the late actinides where 
they are more localized and hence strongly corre- 
lated. However there is currently no consensus on 
the underlying electronic structure of Pu and the 
precise mechanism underlying the anomalous behav- 
ior of plutonium and the late actinides, and several 
efforts in developing realistic approaches to this class 
of problems are currently being pursued. 

The traditional electronic structure techniques 
have been unable to account for the unusual proper- 
ties of Pu. If one supposes itinerant nature of the Pu 
5/ electrons, volume of both the low-temperature a 
phase and that of the fee high temperature 5 phase 
can be described reasonably well [2||. However, a sig- 
nificant spin polarization then results in the LSDA 
calculations [3, 0, H, d, 0]i while a number experi- 
mental investigations have led to conclusion that Pu 
is non- magnetic (for review see Q). 

The strongly localized limit for the 5/ shell has 



been studied by means of the LDA+U technique 
0, 03, EE E3- Depending on the type of the 
" double-counting" correction employed and depend- 
ing on the relative strength of the multiplets and 
the spin orbit coupling and the strength of the local 
Coulomb interaction between 5/ electrons one ob- 
tains either strongly magnetic f 5 or a non-magnetic 
/ 6 configurations for the localized / shell [HI Hit - 
In the latter case one obtains the equilibrium vol- 
ume of the 5-phase in good agreement with the ex- 
periment That has led to conclusion that the 
ground state of <5-Pu should be essentially in the / 6 
non-magnetic configuration, with the 5/ shell being 
split due to the spin-orbit coupling into the com- 
pletely filled J 5 / 2 and empty f 7 ^ 2 states, with the 
local Coulomb interaction U further increasing this 
splitting. Therefore one predicts the ground state of 
S-Pu to be fairly similar to that of Am, though with 
substantially smaller degree of localization of the / 
electrons. Similar results have been obtained re- 
cently by Shick et. al. [3 by means of LDA+DMFT 
using Hubbard I as the impurity solver. 

However, in this picture one expects to have al- 
most no /-electrons in vicinity of the Fermi level, in 
sharp contradiction to the experimental PES [14| . 
where the high peak attributed to / states is clearly 
observed at the Fermi level. That deficiency has 
been addressed in Ref. [l5|, where the dynami- 
cal fluctuations around / 6 ground state have been 
added by means of the dynamical mean-field the- 
ory (DMFT) [16] in conjunction with the pertur- 
bative treatment of the local quantum fluctuations 
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by the fluctuating exc hang e and T-matrix technique 
(T-matrix FLEX) [Hdldi]. The dynamical fluc- 
tuations have not changed the / shell occupancy, 
however they led to essential modification of the den- 
sity of states (DOS), producing a peak which now is 
closer to ( but still separated by a finite energy from 
) the Fermi level, substantially improving the overall 
agreement between the theoretical DOS and exper- 
imental PES for <5-Pu. 

Another school of thought, dating back to Jo- 
hansson [20j suggests that in the solid, plutonium 
is close to an 5/ 5 configuration. These ideas were 
supported by the first LDA+DMFT computations 
in this material by by Savrasov et al. [2l| . At the 
reasonable value of the Coulomb interaction U =4 
eV Savrasov et al. have obtained a double-minima 
shape of the total energy versus volume curve, with 
positions of the minima corresponding well to the 
experimental equilibrium volumes of the a and 6 
phases of Pu, respectively. The minima have been 
assigned to the itinerant and localized state of the 
5/ shell, thus the a — > S phase transition is due to 
the Mott localization of the Pu / electrons. The 
same explanation for the a — > 6 phase transition 
has been put forward by Katsnelson et al. [22| . 
The physical picture underlying this approach, is 
that the actinides provide a complex generalization 
of the localization-delocalization DMFT phase dia- 
gram obtained in simple model Hamiltonians [23| . 
obtained by incorporating the realistic band struc- 
ture, the realistic electronic structure of the f elec- 
trons, and most important, by relaxing the position 
of the atomic coordinates that were kept fixed in the 
model Hamiltonian approach. The LDA+DMFT 
DOS reported in Ref. [2l| did not include realistic 
multiplets structure. Its recent incorporation has 
lead to rather good agreement with the experimen- 
tal PES, and in addition has shown that the fully 
self consistent DMFT solution is non-magnetic [24| . 
by allowing the possibility of ordered states. The 
ground state configuration of both the a and <5-Pu 
were found to be close to the "open-shell" f 5 con- 
figuration, in the sense that the f occupation in the 
ground state is slightly bigger than rif — 5. 

At this point we have, within the Electronic Struc- 
ture Method +DMFT framework, two rather differ- 
ent pictures of the ground state of Pu. In the " near 
5/ 6 " of Refs. one has essentially the 5/ 6 

"closed-shell" non-magnetic configuration, (result- 
ing from a strongly interacting one body Hamilto- 
nian obtained in LDA+U) with dynamical fluctua- 
tions reducing the f occupancy by inducing virtual 
fluctuations into the 5/ 5 manifold, ( and all other f 
occupancies). Another approach of Ref. [2l| can be 
described as starting from an f 5 "open-shell" config- 
uration, which is then screened by the spd conduc- 
tion electrons, what results in a Kondo resonance 
being formed at the Fermi level. The one electron 



Hamiltonian in this case is obtained from the non- 
magnetic weakly correlated LDA. These two phys- 
ical pictures are not orthogonal. In a solid the oc- 
cupation is non integer, and the theory of the An- 
derson impurity model teaches us that at non in- 
teger occupancy the ground state of the model can 
be adiabatically continued from both the / 5 and / 6 
limits. Calculations based on both physical pictures 
produce a ground state and density of states of S-Pu 
in reasonable agreement to experiment. Still, these 
two pictures are substantially different in content, 
and in order to chose which of two approaches (and 
which one of the physical pictures) provides the best 
description of the actual plutonium it is necessary to 
further compare their predictions for other physical 
quantities sensitive to the electronic structure. The 
electronic specific heat coefficient 7 is a convenient 
choice as it is directly linked to strength of the elec- 
tronic correlations. Moreover, it has been measured 
in both a and 5-Pu, and the value of 7 in the 5 phase 
(43 and 64 mJ K -2 mol -1 in S-Pu, stabilized by Ga 
and Al, respectively [1, [25| ; in the range of 35-55 mJ 
K+ 2 mol -1 in the Pug2Amg alloy [26|]) is substan- 
tially larger than in a-Pu (17 mJ K -2 mol -1 [Hj]). 
The specific heat is a measure of the number of de- 
grees of freedom available to the system at a given 
temperature. It is much larger in the case of an 
open shell configuration screened by spd electrons, 
because the degrees of freedom of the / configura- 
tion are transferred to the low energy quasiparticle 
of the system. Hence, in the present work we com- 
pare and analyze predictions of the two picture ( 
/ 5 + and / 6 - ) implemented within LDA+DMFT 
regarding values the 7 coefficient in the a and <5-Pu. 



II. COMPUTATIONAL METHOD 

We have employed the LDA+DMFT technique 
on the basis of the full-potential linear MT-orbitals 
method (FPLMTO) [23] and in conjunction with 
the spin-orbital T-matrix FLEX quantum impurity 
solver [lj|. The details of our technique are de- 
scribed in Ref. We start with the LDA or 
LDA+U calculations of a and <5-Pu by means of 
the FPLMTO method. The calculations of the fee 
5 — Pu phase were carried out at experimental lattice 
parameter (a =4.64 A, V—25 A 3 /at). The complex 
structure of the a phase has been mimicked by an or- 
thorhombicaly distorted diamond structure (which 
has been proved to have almost the same DOS and 
total energy as the actual a-Pu structure [28|) at 
experimental atomic volume (a =3.63 A, 6/a=1.61, 
c/a=2.09, U=20 A 3 /at). We have employed 242 
and 150 k— points in the irreducible Brillouin zone 
for the cubic and orthorhombic structures, respec- 
tively. In both the LDA+U and LDA+DMFT we 
have taken the values of 3 and 0.55 eV for the pa- 
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rameters U and J of the local Coulomb interaction, 
respectively. In the LDA+U calculations we have 
employed the " around mean- field" expression for the 
double counting term (loj . 

The one-particle Hamiltonian H t (k) was obtained 
by Lowdin orthogonalization (using the square 
root of the overlap matrix) of the converged LDA 
or LDA+U calculations. Below we designate as 
LDA+FLEX or LDA+U+FLEX the cases where the 
LDA or LDA+U Hamiltonian, respectively, has been 
employed in the LDA+DMFT calculations in con- 
junction with the spin-orbit T-matrix FLEX quan- 
tum impurity solver. First, the local Green's func- 
tion is obtained by means of the Brillouin zone (BZ) 
integration 

G(iu) = + M)l " H t (k) - ^{itu)]- 1 (1) 

k 

where u> — (2n + l)nT are the fermionic Mat- 
subara frequencies for a given temperature T, fi 
is the chemical potential and T, dc (iuj) is the lo- 
cal self-energy with a "double counting" term sub- 
tracted, viewed as a functional of the local Greens 
function G{iu>) which is solved for self consistently 

Then within the DMFT scheme the local self- 
energy is obtained by the solution of the many-body 
problem for a single quantum impurity coupled to 
an effective electronic "bath" through the Weiss field 
function [16[ 

g- 1 {iLj) = G ff - 1 (iu J ) + Z dc (iuj). (2) 

Gff denotes the ff block of G in equation |T]) In the 
spin-orbit T-matrix FLEX quantum impurity solver 
the local self-energy is obtained as a sum of three 
contributions: 

E = S(™)+S< Tf » + E(™), (3) 

where £( Tff ) and X( TF ) are the T-matrix "Hartree" 
and "Fock" contributions, respectively, T,( pH > is the 
particle-hole contribution. £( Tff ' and £( TF ) are ob- 
tained by substitution in the corresponding Hartree 
and Fock diagrams the bare Coulomb interaction 
with the frequency dependent T-matrix, the latter 
is given by summation of the ladder diagrams in the 
particle-particle channel. E( PH ) is obtained by the 
RPA-type summation in the particle-hole channel 
with the bare vertex being substituted by static limit 
of the T-matrix [II, 0J]. 

For the double counting term in the DMFT com- 
putations we have used the static limit of the self 
energy [l8|, [2^] for the LDA Hamiltonian and the 
Hartree- Fork term for the LDA+U Hamiltonian . 
The DMFT calculations have been carried out using 
the 512-points Matsubara frequencies corresponding 
to temperature 950 K. We carried out DMTF itera- 
tions until convergence in both the chemical poten- 
tial jj, and the local self-energy was achieved. Then 
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FIG. 1: The total (solid line) and partial /-states 
(dashed line) densities of states of a (left panel) and 
5 (right panel) plutonium in the LDA, LDA+FLEX, 
LDA+U and LDA+U+FLEX approaches, respectively. 



the Pade approximant was used [301 ] for analytical 
continuation of the local self-energy to the real axis 
in order to obtain the DOS n(E) = -1/ttTiQ[G{E)]. 
The electronic specific heat coefficient 7 was com- 
puted as ^-k%tr[N(E F )Z- 1 ], where N{E F ) is the 
DOS matrix at the Fermi level (which for rather 
low temperature used here can be approximated 
to /i) and Z^ 1 is the quasiparticle residue matrix, 
Z- 1 = I - d$[E(E)}/dE\ E=0 . where I is the iden- 
tity matrix, and Sigma is zero outside the f block. 
The average mass enhancement ^ is defined by 

tT ^tr[Nf]^ ' wnere Nf is the f spectral function ma- 
trix at the Fermi level. 

III. RESULTS 

The calculated total and partial /-states LDA and 
LDA+U DOS are presented in Fig. [I] Already in the 
case of LDA one may clearly distinct the filled / 5 / 2 
and empty / 7 / 2 states on the DOS of S-Pu. For a- 
Pu the separation between the / 5//2 and f 7 ^ 2 states 
is somewhat smeared out, as has been also observed 
in Ref. [TTJ] . However, in both the a and S phases 
there is still substantial density of the /-states at 





Effective mass 
a-Pu <5-Pu 


7 (mJ KTVo-r 1 ) 
q-Pu <5-Pu 


LDA+FLEX 
LDA+U+FLEX 
Exper. [8, 25, 26] 


3.85 4.02 
2.22 1.68 


20.7 37.9 
15.6 8.6 
17 43, 64, 35-55 



TABLE I: Effective mass — and the electronic specific 
heat coefficient 7 



the Fermi level Nf(Ep), therefore within LDA one 
obtains an "open-shell" configuration of the /-band. 
The LDA+U calculations of 5-Pu lead to familiar 
picture [Tol [TP . [ToT ] of non-magnetic configuration 
with the filled / 5 / 2 and empty f 7 ^ 2 states, well sep- 
arated in energy by value of U, negligible Nf(Ep) 
and the / shell occupancy close to 6. In a-Pu sub- 
stantial broadening and intermixing of the / 5 / 2 and 
J 7 / 2 states persist in the LDA+U, retaining appre- 
ciable N f (E F ). 

The values of the electronic specific heat coef- 
ficient 7 and effective mass m* calculated within 
the LDA+FLEX and LDA+U+FLEX techniques 
are listed in Table HI One may notice, that mod- 
erate values of the effective mass m* in the case 
of LDA+U+FLEX indicate rather weak correlations 
present in the system. In fact, the effective mass m* 
in S-Pu appears to be somewhat smaller than in a- 
Pu. Similarly, in LDA+U+FLEX the 7 coefficient 
of specific heat in the S phase is nearly two times 
smaller as compared to a-Pu, contrary to experi- 
mental measurements, which shows 7 bigger than 
two times and sometimes almost three times larger 
in 5-Pu. Almost order of magnitude difference be- 
tween experimental and theoretical 7 in the case of 
S-Pu means that in the "closed-shell" configuration 
of 5-Pu electronic correlations appear to be much 
weaker than one would expect on basis of exper- 
imental evidences. This should be contrasted with 
the LDA+FLEX calculations which predict effective 
mass in 5-Pu to be quite large (4.02) and somewhat 
larger, than in a-Pu. The LDA+FLEX technique 
also reproduce well experimental tendency of en- 
hancement of the value of 7 in the 5 phase as com- 
pared to a-Pu. Actual values of 7 obtained by means 
of the LDA+FLEX technique are rather close to ex- 
perimental results in both the a and 7 phases. 

In Fig. [2] we display the self-energy T, mm (iu)) (for 
m—-i) computed on the Matsubara axis within the 
LDA+FLEX and LDA+U+FLEX techniques for a 
and <5-Pu. One may notice that 3E(iu;) obtained 
by LDA+U+FLEX is about four times smaller than 
one calculated within LDA+FLEX. That should be 
expected on the basis of already observed tendency 
for m* . We stress again, that in both descrip- 
tions (LDA+FLEX and LDA+U+FLEX ) we have 
a quasiparticle feature in the spectra. This feature, 
becomes the famous Kondo resonance as U increases. 
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FIG. 2: The imaginary part of self-energy E mm (iw) for 
m=-3 on the Matsubara axis, computed for a and 5- 
Pu within the LDA+FLEX and LDA+U+FLEX tech- 
niques. 
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FIG. 3: The theoretical LDA+FLEX and experimental 
PES [3 for a and <5-Pu 



However, the values of the many body renormaliza- 
tions , namely the value of the slope of SS(icj) ob- 
tained in the LDA+U+FLEX calculations of <5-Pu is 
relatively small and closer to the non interacting sys- 
tem (see Fig. [2]) while the peak at Ep that is present 
in the LDA+FLEX DOS is closer to the Fermi level 
and more enhanced. 

Finally in Fig. [3]we compare the a and 5-Pu DOS 
obtained within the "open-shell" LDA+FLEX ap- 
proach to the experimental photoemission spectra, 
collected at a photon energy of 40.8 eV, where the 
band 6d and 5/ cross section are nearly equal [T3 |. 
In order to produce the theoretical spectra we have 
multiplied the total DOS by the Fermi function, cor- 
responding to experimental temperature of 80 K and 
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normalized it to equal spectral weight between -4 
and eV for a and <5-Pu as was done in analysis of 
the experimental PES in Ref. [3| ■ In result one has 
very good agreement between experimental and the- 
oretical PES in vicinity of the Fermi level. For lower 
energies below -0.5 eV the agreement is less satisfac- 
tory, but there the d states should give appreciated 
contribution to PES, therefore the accurate evalua- 
tion of the corresponding 6c? and 5/ matrix elements 
may become necessary for a faithful representation 
of PES. Moreover, the multiplet effects in the 5/ 
shell should be taken into account and these are be- 
yond our simple T-matrix FLEX impurity solver. 

In conclusion, we have carried out comparative 
analysis of the specific heat and electronic structure 
in a and <5-plutonium by means of the LDA+DMFT 
technique in conjunction with the perturbative T- 
matrix and fluctuating exchange quantum impu- 
rity solver. We have shown that by assuming the 
"closed" / 6 configuration of the Pu / shell one can 
not explain high value of the electronic specific heat 
coefficient 7 in S-Pvl as well as substantial enhance- 
ment of 7 in the S phase as compared to a-Pu, 
even in the case when the dynamical fluctuations 
are properly taken into account. At the same time, 
by assuming rather the "open shell" configuration, 
which is closer to / 5 one may obtain the observed 
enhancement of 7 in the 8 phase and actual values 
of 7 in rather good agreement with the experiment. 
The experimental PES is also well reproduced by 
the LDA+DMFT DOS in the "open shell" configu- 
ration. This suggests that the actual Q ground state 



of the / shell in Pu should be rather closer to the 
"open shell" / 5 configuration than / 6 . This picture 
is consistent with the results of a series of EELS and 
XAS studies of Ref. [3l[ and in goo d ag reement with 
the results of recent LDA+DMFT [24[ calculations, 
where more sophisticated solvers were used to ac- 
cess the strong correlation regime. The good agree- 
ment with respect to experiments achieved by our 
work suggest the possibility of exploring the param- 
agnetic state of inhomogeneous alloys and interfaces 
of actinides using DMFT LDA+FLEX. Notice in the 
context that simplified solvers such as the iterated 
perturbation theory was very useful for the study of 
the paramagnetic phase of the one band Hubbard 
model, but far worse in the description of magnet- 
ically ordered phases. For this reason we have not 
addressed in this paper the issue of magnetic long 
range order. This issue is very delicate, because as 
is well known in the context of the Anderson im- 
purity model, approximate perturbative treatments 
that signal the onset of magnetism, can indicate the 
need to include explicitly Kondo screening rather 
than the occurrence of a magnetic instability. 
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